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We describe the gluon plasma in zeroth order as a gas of free quasi-particles with a temperature- 
independent dispersion relation of Gribov type, E (k) = \J k 1 + ^£ , that results from the restriction 
of the physical state space to the Gribov region. The effective mass controls infrared diver- 
gences of finite-temperature perturbation theory. The equation of state of this gas is calculated 
and compared with numerical lattice data. At high temperature T there are power corrections to 
the Stefan-Boltzmann law that are of relative order 1 /T 3 for pressure and 1/T 4 for energy. 
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1. Introduction 

The equation of state (EOS) of the gluon plasma is rather precisely known from numerical 
studies in lattice gauge theory [p, and there exist excellent phenomenological fits to the EOS 
of the gluon plasma in the deconfined phase In contrast with this practical success, it was 
discovered by Linde [^] as far back as 1980 that finite-temperature perturbative QCD suffers at a 
fundamental level from infrared divergences, which suggests that finite-temperature perturbation 
theory neglects an essential feature of QCD. It has been proposed [Q] to control these divergences 
by introducing a magnetic mass m ~ g 2 T. It has also been proposed ^ that the divergences stem 
from an inadequate application of the principle of gauge equivalence. Indeed in 1978, two years 
before Linde 's discovery of the infrared divergence [B], Gribov showed that infrared modes are 
strongly suppressed when gauge equivalence is imposed at the non-perturbative level [Q]. 

Two or more different configurations may be gauge equivalent even though both satisfy the 
linear Coulomb gauge condition, £^ =1 9,A,- = 0. When enumerating physical states, only one of 
these "Gribov" or gauge copies should be counted, so the space of physical states is reduced to the 
fundamental modular region (FMR), a region that is free of Gribov copies. Gribov [Q] found that 
the dispersion relation E(k) = k gets modified because of the reduction of the physical state space, 
and he obtained instead 

E(£) = ^2 + _ (L1) 

where k = |k| , and M is a QCD mass scale. The reduction of the physical state space was originally 
proposed as an essential feature of the confinement mechanism [|], [7|]. However statistical mechan- 
ics is primarily a matter of counting states, and the reduction of the physical state space required 
by the gauge principle influences the EOS at all temperatures. Here we shall be concerned with its 
effect in the deconfined phase. 

It is known from numerical studies ^ that at high temperature the EOS of the gluon plasma 
approaches the Stefan-Boltzmann law, £ = 3p = 3c$bT 4 , where £ is the energy per unit volume, p 
is the pressure, T is the temperature, and csb = 7jfi{N 2 — 1) in SU(N) gauge theory. Thus it seems 
reasonable to describe the gluon plasma at high temperature in first approximation as a gas of non- 
interacting quasi-particles. We shall describe the quasi-particles by the Gribov dispersion relation 



( |l.l| ) or a similar one, for E(k) is only approximately known. Fortunately the results obtained hold 
under rather general conditions on E(k). We call the gas of non-interacting quasi-particles, with 
modified dispersion relation, the FMR gas. 

Many aspects of the Gribov scenario have been verified in recent investigations. Infrared 
suppression of the gluon propagator in Coulomb gauge has been observed in numerical simulation 



], but less strongly in [W, and in Landau gauge in 3-dimensions [10]. It has also been found in 



variational calculations in Coulomb gauge Ql 1|], and in Schwinger-Dyson calculations in Coulomb 
[12] and in Landau [13] gauge. A long-range color-Coulomb potential was found in numerical 
simulations in Coulomb gauge in the deconfined phase [|l4|], and is reported in the talk by Stefan 
Olejmk at this conference. 

The questions we wish to address here are: (i) What is the equation of state of the FMR gas? 
(ii) How does it compare with the EOS that is known from numerical studies? (More details are 



provided in ||15|].) 
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2. EOS of the FMR gas 

In the Stefan-Boltzmann limit, the degrees of freedom for each gluon momentum k are the 
two states of polarization, each with color multiplicity (N 2 — 1). These are precisely the degrees 
of freedom in Coulomb gauge, and we shall use this gauge for our calculation. It is a "physical" 
gauge without negative metric states, and all constraints are satisfied identically. Although the 
Coulomb gauge is not manifestly Lorentz covariant, this is not necessarily a disadvantage in the 
deconfined phase because at finite T the heat bath provides a preferred Lorentz frame, and the 
manifest symmetries of the Coulomb gauge are the symmetries of the physical problem at hand. 

The partition function of the FMR gas is given by the Planck distribution, Z = \\ n [ 1 — 
exp(— $E n ) where [3 = l/T is the inverse temperature, and n = (k,X,a), where k is 3-momentum, 
X = 1,2 is polarization, and a = l,...(N 2 — 1) is color. For the energy density, £ = one 
obtains 

(N 2 -l) r k 2 

L > I dk --— E(k), (2.1) 



7i 2 Jo exp[p£(ifc)] - 1 
and for the pressure, p = T^^- = ylnZ, 



37i 2 Jo exp[$E{k)} - 1 dk 
and for the trace anomaly 6 = £ — 3p, 

(N 2 -l) r°° d ( -E{k) \ 

% 2 Jo exp[p£(*)]-i a*v * r K ] 

We also have £ = T 2 -^{^), from which it follows that the trace anomaly may be written = 
£ — 3p = T 5 -^^). Upon integration this yields 

ooo 

p = CsB T 4 -T 4 dT' T'~ 5 e(r'). (2.4) 



T 

3. FMR gas at high temperature 

Suppose that the leading deviation of E(k) from k at high k is expressed by a power law, 

E(k)/k = \+c/tf+... . (3.1) 



For the Gribov dispersion relation E(k) = v / k 2 +M 4 /k 2 one has y = 4, whereas for a gluon mass, 
E (k) = Vk 2 + m 2 , one has y = 2. The gluon condensate has dimension 4, which leads one to expect 
y = 4, whereas if there were a condensate of dimension 2, one would expect y = 2. The asymptotic 
behavior of the EOS is qualitatively different for y greater or less than 3. We suppose y > 3, and 
consequently the deviation from the ultraviolet behavior E(k) = k is soft. 

For y > 3, the asymptotic high-J limit of the trace anomaly is obtained from the substitution 



exp[^^] — 1 — > m fl23| ). This substitution cannot be made in the integrals for p and £ because 



they would diverge. It gives a linear asymptotic trace anomaly, 



8 = 17 + 0(1), 



(3.2) 
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where 



(AT 2 -1) f k 4 d (~E(k) 



7i 2 Jo E(k) dk\ k 

p CO 

= 3(N 2 - 1)k- 2 / dkk 2 \n[E(k)/k], (3.3) 
Jo 

is an integral that converges for y > 3 (as we have supposed), and is positive for E(k) > k. These 



are sufficient conditions for the linear asymptotic form ( p.2| ) with L > 0. Although L describes the 
high-r limit of the FMR gas, the last integral involves all momenta k. For the special case of the 
Gribov dispersion relation one obtains L = (N 2 — 1) {%\fl)~ l M 3 , which is proportional to 
M 3 although only M 4 appears in the dispersion relation. 



The pressure at high T is obtained from ( |2.4j ) which yields 

p = c SB T 4 -(l/3)LT + 0(l), (3.4) 

Thus for the FMR gas, the leading deviation of the pressure from the Stefan-Boltzmann law is 
linear in T. However this linear term — and only a linear term — is annihilated in the formula for 
the energy density £ = r 2 ^(^), which gives 

e = 3c SB T 4 + 0{l). (3.5) 

An EOS of this type was obtained as a fit to the lattice data in [ |l6| ] . 

4. Comparison with numerical EOS 



The EOS of the FMR gas at high T is not sensitive to the exact form of E(k) because (3.4) 



and (3^5) hold as long as y > 3 and E(k) > k. To compare with the numerical data, we take the 
Gribov dispersion relation The unknown mass scale M is determined by fitting the anomaly 
6 = £ — 3p at high T, because the corrections to this quantity are expected to be small. [Indeed, 
taking thermal perturbation theory as a guide [|17|], we note that the leading correction to p is 
of order Ap ~ g 2 (T)T 4 ~ ^y. The anomaly is given by 9 = r 5 ^(^), so the corresponding 



correction to the anomaly, A6 = r 5 ^(f£), is of order T 5 ^( ^ = ~ g 4 (T)T 4 .] 

The numerical data of ^ for ^ are represented by the black interpolating curves in Fig. 1 (left). 
(More recent studies include dynamical fermions or a chemical potential that cannot be described 
by the FMR gas.) The data for N x = 6 and 8 agree, and were interpreted as continuum values [j|]. 



The red dots are obtained from the analytic formula (2.3), with mass scale set at M = 2.6T C by 
fitting at high T, where T c is the transition temperature. The relatively large difference in the tran- 
sition region between the FMR gas and the numerical data for 9 = £ — 3p occurs because the FMR 
gas does not exhibit a sharp phase transition, whereas for pure SU(3) gauge theory there is a first 
order phase change, so £ is discontinuous while p is continuous. We do not attempt to estimate the 
error of M because perturbative-type corrections to the FMR gas have been neglected. From the 
value T c = 0.625-y/a of [JIJ], where a is the string tension, one gets M = l^y^, or M = 705 MeV, 
where the string tension for the quarkless theory is defined to be ^/o = 440 MeV. 

The numerical data of [|l|] for jj, ^ and |§ are are displayed as black interpolating curves in 
Fig. 1 (right), where jj = ^f. The horizontal line represents the Stefan-Boltzmann EOS. For the 
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Figure 1: Numerical and analytic plots of (e- 3p)/T 4 (left), and of e/T 4 , 3/4 s/T 3 , 3p/T 4 (right). 



FMR gas, |§ and ^ approach the Stefan-Boltzmann limit like jj and respectively, being quite 
close to it at the highest temperature displayed, T = 5T C , whereas the gluon plasma approaches the 
Stefan-Boltzmann limit more slowly. The difference between the FMR gas and the gluon plasma 
in the range 2T C to 5T C appears attributable to perturbative-type corrections of moderate size. In 
standard thermal perturbation theory [17] these are of leading order g 2 (T) ~ j-^, but diverge at 
order g 6 , whereas corrections to the FMR gas are expected to be calculable. 



5. Discussion 



Above the transition region, the EOS of the FMR gas gives a good description of the most 
prominent feature of the gluon plasma which is the rapid drop of the pressure compared to the 
energy from the Stefan-Boltzmann value, as T decreases from infinity. The linear asymptotic trace 
anomaly ( |3~^ ) provides a ready explanation for this, that holds for any quasi-particle model with 



y > 3 and y^jp-j < 0, although other fits are certainly not excluded. The FMR gas is not exact 
even at high temperature because of perturbative-type corrections. We expect however that they 
are calculable and of moderate size above the transition region. 

The transition region is not so well described by the FMR gas. There is no shatp phase tran- 
sition because p(T) is an analytic function. Moreover the dependence on N [of SU(N)] is only 
through the coefficient (N 2 — 1), whereas even the order of the phase change depends on N, being 
second order for SU(2) and first order for SU(3). An analysis of the phase transition based on cen- 



ter symmetry is given in [18]. However the dispersion relation E(k) = y k 2 + ^ has a minimum 
energy E m j n = v2M, so the thermodynamic functions £, p and s of the FMR gas are exponentially 
small for T < V2M = 997 MeV (for M = 705 MeV). The mass of the lightest glueball is of order 
1 GeV, so the thermodynamic functions of the FMR gas are exponentially small where they are 
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supposed to be. Altogether the FMR gas, with a single parameter M which is the mass scale, is 
competitive with cut-off phenomenological models 

The improved phenomenological models of [Q] are more precise than the FMR gas especially 
in the transition region, but they also require several parameters, whereas we have fit only the QCD 
mass scale M. However the FMR gas is not intended to be a precise phenomenological model, but 
rather to provide a useful starting point, well founded in the principles of gauge theory, that allows 
calculable, moderate size corrections at high T . Athough it is denned by Gribov's dispersion rela- 
tion of 1978, the FMR gas has two important properties that were later independently reinvented: 
(i) Its EOS closely resembles simple phenomenological models of the gluon plasma [|l^]. (ii) The 
effective mass ^- controls infrared divergences in higher order corrections so it is not necessary to 
introduce the magnetic mass m ~ g 2 T [Q] for this purpose. 
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